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Abstract 



The nonmesonic decay of the hypertriton is calculated based on a hyper- 
triton wavefunction and 3N scattering states, which are rigorous solutions 
of 3-body Faddeev equations using realistic NN and hyperon-nucleon inter- 
actions. The pion-exchange together with heavier meson exchanges for the 
AA^ NN transition is considered. The total nonmesonic decay rate is found 
to be 0.5% of the free A decay rate. Integrated as well as differential decay 
rates are given. The p- and n- induced decays are discussed thoroughly and it 
is shown that the corresponding total rates cannot be measured individually. 
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I. INTRODUCTION 



The hypertriton consisting of a neutron, proton and a A or S, which strongly convert into 
each other, is bound against A-deuteron decay by 0.13 ± 0.05 MeV. Recently that number 
could be reproduced [Q] by solving the Faddeev equations based on realistic NN forces and 
the Nijmegen hyperon-nucleon interaction [Q. The resulting wavefunction has all two-body 
correlations exactly built in as enforced by the various two-body forces. As the lightest 
hypernucleus the hypertriton plays the same role in hypernuclear physics that the deuteron 
does in nuclear physics. However, in contrast to the deuteron, the hypertriton is subject 
to the weak decay and has a lifetime comparable to that of the free A, ta = 2.63 x 10^^° 
seconds. The first data on light hj^ernuclear lifetimes have been obtained using bubble 
chamber experiments and emulsion works which in most cases only detected the mesonic 
decay modes. These measurements suffered from low precision, very poor statistics and 
difficulties with particle identification, leading to fairly large error bars. Along with the 
mesonic two-body decay mode \H {71^)+^ He{^H) there are the corresponding mesonic 
multi-body decay modes \H 7r^(7r*^) + d + p{n) and \H 7i^{7i^) + p + n + p{n). The 
most precise experiment to date for the combined two- and multi-body decay modes gave a 
value of r = (2.28 -|- 0.46 — 0.33) x 10~^° seconds 0. A similar measurement furthermore 
was also able to estimate the decay branching ratio T{\H — > 7i~ +^ He)/T{\H — > all 7i~ 
mesonic modes) as 0.30 ± 0.07. 

Besides the mesonic decay channels there are also two nonmesonic modes, \H d + p 
and \H p + p + n. While it is well known that these channels dominate the weak decay of 
heavy hypernuclei they are expected to be rare for the hypertriton since the mesonic modes 
are not Pauli suppressed [Q. As a first step this two-baryon decay mode AA^ — > A^A^ can be 
understood in terms of the free-space decay mechanism A nN with virtual pion that is 
absorbed on a second nucleon bound in the hypernucleus . However, the large momentum 
transfer involved in the reaction leads to a mechanism that is sensitive to the short distance 
behavior of the amplitude and allows for the exchange of heavier mesons. The production of 
these mesons would be below threshold for the free-space A decay, but they can contribute 
through virtual exchange in a two-baryon decay channel. 

The weak nonmesonic decay channel is important since it allows access to the fundamen- 
tal aspects of the four-fermion, strangeness changing weak interaction. In general, starting 
with the Standard Model electroweak Hamiltonian and taking into account QCD corrections 
at short distances yields an effective weak V — A interaction that could presumably pre- 
dict the relative strength of the A5'=0 and AS=1 transition. Thus, hadronic weak matrix 
elements of the form {MB'\Hw \ B) can be calculated [^]. Using these weak vertices as 
a starting point for effective nuclear two-body operators that are then implanted into the 
nucleus with the usual nuclear many-body wave functions provides the testing ground for 
the effective interaction. 

Parity violation in hadronic systems represents a unique tool to study aspects of the 
nonleptonic weak interaction between hadrons. The nonmesonic process resembles the weak 
AS=0 nucleon-nucleon interaction that has been explored experimentally in parity-violating 
NN scattering measurements. However, the AA^ NN two-body decay mode can explore 
both the parity-conserving (PC) and the parity-violating (PV) sector of the AS=1 weak 
baryon-baryon interaction while in the weak NN system one is limited to the weak PV 
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interaction. A number of theoretical approaches to the AA^ NN decay mode have been 
developed over the last thirty years which are more extensively reviewed in Ref. f^. The 
AS" = weak nucleon-nucleon interaction at low and intermediate energies has generally 
been described in a meson exchange model involving one strong interaction vertex and 
one weak one; the same approach can be used for a microscopic description of the AS=1 
AN NN mechanism. 

A recently completed major study of the nonmesonic decay of p-shell hypernuclei 
found that proper short-range correlations in the initial and final state are of major im- 
portance in predicting decay rates and asymmetry observables. However, in a shell-model 
framework bound state wave functions, spectroscopic factors, short-range correlations and 
final state interactions do not all originate from the same underlying dynamics and, there- 
fore, introduce approximations that may be difficult to quantify. Since the aim of investi- 
gating the nonmesonic decay is to extract information on hadronic weak vertices from the 
AA^ — i> A^A^ process the decay of few-body hypernuclei offers a window since all nuclear 
structure ingredients are derived from the same baryon-baryon interaction. 

It thus appears worthwhile to repeat a former study 0] on the weak nonmesonic decay of 
the hypertriton, where a simplified uncorrelated deuteron-A wavefunction has been used. We 
expect that correlations should play an important role, since the mesons emitted by the weak 
hyperon-nucleon transition are reabsorbed by the nucleons. The resulting meson-exchange 
operator acts like a two-body force and consequently probes the hypertriton wavefunction 
in its dependence on the pair distance between a hyperon and a nucleon. Furthermore, the 
final three nucleons will interact strongly with each other, which might infiuence significantly 
the decay process. This dynamical ingredient has been neglected in Ref. and will now be 
fully incorporated. 

In section |T| the theoretical formalism for the evaluation of the decay matrix element will 



be given. Section [II describes the necessary technicalities. A special section IV is devoted 



to the meson-exchange operator. We present our results in section [V]. We summarize and 



conclude in section VI 




II. FORMALISM 

There are two nonmesonic decay channels 

According to standard rules, the partial decay probabilities in the total momentum zero 
frame are 

= - y \< I o I ^3^„ >p 

2 ^ ' kNkarriNma ' ' ' 



2Mn 2Md 
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and 

2 ^-^ ' ki k2 ks mi 1712 ma ' 
m mi m2 ms 



dk, a, 2. i + + 4,) i I M. „ - 3M. - ^ - ^ - ^ I , (2) 

where are appropriate three-nucleon scattering states, O the transition operator and 
\l/3 fj the hypertriton wavefunction. We use throughout nonrelativistic notation. The binding 
energies are defined as usual in terms of various masses as 



Ms^H = SMjv + Ma + e 

Ma = 2Mn + ea (3) 

Further we introduce Jacobi momenta for the final 3N states 

— * — * 

q=l{h-\ ih + 4)) (4) 

— * — * 

and identify for the nd breakup — and ki + k2 — k^- Then some simple algebra leads 
to 

dr^'^' = 1 E \< "^tLm, I O I ^3^^ >r (5) 
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with 



go = y ^(Ma - + e - ed). (6) 

Because of the averaging over spin directions the matrix element squared is independent of 
^0 and we get just a number 

"-' mmj^m^ 



Similar steps lead to 

rfr"+"^+^ =\ I< 4;lim.m3 I O I ^iHm >? 271^ qdqdpp'dp (8) 

m mi m2 ma 



9 z ^ I pgmim2m3 i i a ' ' — Q 

with 
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Again the spin-averaged part depends only on the angle 9 between p and g, thus 
dr"+-+P = IGTr^^M^gp^rfpsin^rf^i ^ l< ^^Urn.^a I ^ I >P (10) 

^ m mi 1712 1713 

This form is convenient for the integration to determine the total (nnp) decay rate. For the 
display of the angular and energy distribution of the three nucleons, the following, equivalent 
form is more convenient 

m mi 7712 ms 

27rdhdk2dS-^ M^K.K, ^^^^ 

^kf{2k2 + ki ■ k2Y + kl{2k^ + 4 ■ k? 

Here ki and /c2 denote the directions of two detected nucleons. That choice of four angles 
relates the lab energies Ei and E2 kinematically through energy and momentum conserva- 
tion, as shown in the example of an interparticle angle B12 = 180° in Fig. ^ A point on 
that curve can be defined through the arclength S measured from some conveniently chosen 
point. Our choice for S" = is shown in Fig. |l|. Thus instead of expressing the fivefold 
differential cross section with respect to dEi, for instance, we have chosen dS in equation in 
Eq. (p. 

In (0) the necessary additional information, whether the detected particles 1 and 2 are 
a proton-neutron pair or two neutrons, has been dropped for the sake of a simple notation. 
The final scattering state carries, however, additional isotopic spin quantum numbers vi V2 t'a. 

Throughout the paper the normalisation of the momentum states is always like 
<k\k' >= 5^{k - k '). 

The operator O is of two-body character and acts between the A and a nucleon. In the 
hypertriton wavefunction let us denote the A to be particle 1, then 

O = E 0{l,i) (12) 

i=2,3 

Because of the antisymmetry of the hypertriton state and the scattering states with respect 
to the two nucleons 2 and 3 the nuclear matrix element simplifies to 

< I 6 I ^3 ^ > = 2 < I 6(1, 2) I ^3 ^ > (13) 

The exact inclusion of the final state interactions among the 3 final nucleons can be 
performed in analogy to electron scattering on ^He [0]. We exemplify it for the nnp breakup 
process. For our notation in general we refer to [P^. 

The scattering state \E'^"^ = is Faddeev decomposed 

^(-) = {1 + P)iIj(-\ (14) 

where P is the sum of a cyclical and anticyclical permutation of 3 objects and ip^'^ is one 
Faddeev component. It obeys the Faddeev equation 
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(-) 



(15) 



with 



(16) 



P12) \(po>= ^(1 - P12) \p>\q> 



(17) 



Here Gq is the free three-nucleon propagator, t*-"-* the NN (off-shell) t-matrix and takes 
care of the identity of the three nucleons. Note that P12 acts in the two-body subsystem 
described by the relative momentum p. 

Let us now insert (p!^, (|T5|) and (0), into the nuclear matrix element 

< I 6(1, 2) I ^3 > = < I (1 + P)6(l, 2) I > 



< 



:i + P)6(l, 2) I > + < I PtGo(l + P)6(l, 2) I ^3^ > 



< 0S I (1 + tGo){l + P)6(l, 2) I ^3^ > + < 0g I (1 + tG'o)PtGo(l + P)6(l, 2) I ^3^ > 



+ <ro\ il + tGo)PtGoPtGoil + P)dil,2) \^3H> + ■■■ (18) 



In the last equality we iterated (|T5|). It is then easily seen that this can be put into the form 
<vI/(rJ|6(l,2)|vl/3^> = <0SI(l + i^)O(l,2)|vl/3^> + <0g|(l + P)|f/>, (19) 
where | U > obeys the Faddeev equation 

\U> = tGo{l + P)6(l, 2) I ^3 ^ > + tGoP I U > (20) 

Note the driving term of that integral equation contains the operator 6(1, 2) applied to the 
hypertriton bound state | "^j^H > ^ind it also includes rescattering terms of first order in t. 
The pure plane wave impulse approximation is the first term on the right hand side of (P^. 
A similar reduction yields for the nd breakup process 

<4;M0(1,2) |vl/3^> = <0| (1 + P)6(1,2) |M/3^> + <0|P|f/>, (21) 

where now | > contains the deuteron state | yj^ >: 

\(f)> = \^d>\qo> (22) 

and the same state | U > appears. 

Using the weak transition operator 0(1,2) in the context of a strictly nonrelativistic 
framework requires some approximations, which we would like to describe in the example 
of pion exchange . According to Fig. ^ the transition operator is 
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0(1,2) = 75^(^2) f ^^"^o u{ki)GFml {A^ + B^j,) u^{h) (23) 

Here u and -u are the usual Dirac spinors, (^nntt the strong couphng constant for the NNtt 
vertex, and Gpfn'^ = 2.21 x 10^'' the weak couphng constant. The constants A.,^ = 1.05 
and i?7r = —7.15, which determine the strength of the parity violating and parity conserving 
amplitudes, respectively, have been adjusted to reproduce the decay observables of the free 
A particle 0. We assume the same form factor F at the two vertices (the strong and the 
weak one). In the nonrelativistic reduction at the weak vertex the nucleon mass, Mjv, and 
the A mass. Ma, are replaced by their average, M. Then one finds 

0(1. 2) _ -G,rnl |^ ■ i (a, + • ) (24) 

with 

q-ir = h- fc/ = - k2 (25) 

We have to use two types of Jacobi momenta, one set referring to the hypertriton com- 
posed of ANN and another set for the final state of three nucleons. The latter ones have 
been already defined in (^) and will be denoted by p and q . The ones for the hypertriton 
are 

^ _ - MAk2 



-.^ {Mn + Ma)/c3 - MNjki + fca) 

^ 2Mn + Ma ^ ^ 

Then for total momentum zero and using the spectator condition q = q' one has 

'^-=^-^ + 2(Ma + M^) ^ (2^) 

As in the derivation of 0(1, 2) itself we also neglect here the difference M\ — M]y with respect 
to Ma + Mat and put 

q„ — > p-p' (28) 

Then we get an ordinary two-body force, which does not depend on the momentum of the 
third particle (which it would if the mass difference would be included). 

A final remark refers to the isospin part of the transition matrix element. At the weak 
vertex the A has to change into a neutron or a proton by emission of a vr" or vr", respectively. 
This can be formally accomplished by setting artificially the A state to be | | — | > in isospin 
and introducing r at the vertex. This is a well known trick ^,|^ and is in agreement with 
the empirical AI = | rule. As a consequence the two-body force (|2^ ) has to be multiplied 
by fi ■ fa. 
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Now in the hypertriton the A-particle is treated as a strongly interacting particle and 
has therefore isospin zero. The isospin part of the hypertriton (A part only) is 

l^> = l(^^)0>23|00>i, (29) 

where the indices denote the particles and (| |)0 the obvious isospin coupling for the two 
nucleons. 

Now the action of ti ■ T2 resulting from the weak transition requires the A-particle to 
be treated as | ^ ~ | >i consequently the isospin part of the hypertriton has to be 
reinterpreted as 

I ^ > ^1 ^ >wcak = I (2 2)*^ >23| 2 ^ 2 

1,1 1 ,11 ,1 1 1,1 1 ,1 1 ,11 
= >1 >2 >3 7= >1 >2 >3 , (30) 

y/2 2 2 '22 '2 2 ^/2 2 2 '2 2 '22 ^ ^ 

which displays the partner (nucleon 2) of the A to be a proton or a neutron, respectively. 

In the nuclear matrix element one acts from the left by the isospin state of the final three 
nucleon system and the fi ■ T2 operator and gets for the isospin part alone 

1 /3 

3 < {t-)TMT \n-r2\e >„eak = 1 ^ (v^^tO + ^tl) (31) 

The index 3 on the bra indicates that the isospin t refers to the (12) subsystem. 

We see that only total isospin T = ^ occurs for the three final nucleons. For isoscalar 
meson exchanges ti-T2 is replaced by the unit operator and the corresponding matrix element 
is 
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< {t^)TMT I 9 >„eak = i {-^tO + VS^tl) (32) 



One can artificially separate the contributions from the proton and neutron induced decays. 
This corresponds to the first and second parts on the right hand side of (|5(]|), respectively. 
Keeping only the one or the other both isospins T = ^ and T = | contribute, therefore 
(0) and (0) will be adequately modified. That separation into p-induced and n-induced 
decays will be considered in section |V]. 



III. TECHNICALITIES 

The hypertriton state contains the ANN and the T,NN parts. The A — S conversion 
is crucial for the binding of the hypertriton, nevertheless the TiNN admixture is extremely 
small [|l|. Thus we neglect the contribution of the S decay and keep only the ANN part. 

In [U the hypertriton state has been determined in a partial wave representation and 
we refer to for the details of our notation. Here we need only the form 

I >= ^PP^ j d<l<f I PQOt > ^a(M), 
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where p, q are the magnitudes of the Jacobi momenta ( ]26|) and a denotes the following set 
of discrete quantum numbers 

« = (^^)j(A^)/(j7)j(t^)r 

Here {ls)j describe the coupling of orbital angular momentum / and total spin s to the total 
two-body angular momentum j of the AA^ subsystem, (A|)/ the corresponding coupling of 
orbital and spin angular momentum of the other nucleon to its total angular momentum /, 
(jJ) J, the resulting jl coupling to the total angular momentum J and finally the isospin 
coupling of t = I and | to total isospin T = 0, as described above. 



Also for the evaluation of the matrix elements (0) and ( pT]) and the solution of the 
Faddeev equation (|20|) we work in a partial wave representation, using a complete set of 
basis states now for three nucleons. They are again denoted as | pqa >n but adding a 
subscript A^ to indicate that the Jacobi momenta are now from (^). Furthermore one has to 
note that this is a subset of states antisymmetrized in the subsystem of particles 1 and 2, 
thus (/ + s + t) has to be odd. 

Now projecting the Faddeev equation into the basis | pqa >n and inserting appropriate 
decompositions of the unity one gets 

AT < pqa \ U > = 

< pqa I tGo{l + P) I pqa >n n < p'q'a \ 0(1,2) | p"q"a" > + 



t 



N 



< pqa I tGoP I p'q'a' >n n < p'q'a' \ U > (33) 



This is a coupled set of integral equations, with a kernel part, which is well known ^ from 
3N scattering, and an inhomogeneous term, whose part left of 0(1,2) is also familiar from 



electron scattering |]10[. What is left as a new structure is the application of the 0(1,2)- 
matrix onto the wavefunction component of the hypertriton. 

Now that 0(1, 2)-matrix is obviously diagonal in the quantum numbers of the spectator 
nucleon: 

N < pqa \d{l,2)\ p'q'a' >= ^^^^7^ 5a A' 5/ /' < p(/s)j | 6(1, 2) | p'(/V)j > (34) 

qq' 

and one is left with a simple application of the two-body force onto the hypertriton in 
momentum space. The right hand side of ( P^ should contain the appropriate isospin matrix 
element in the 3-particle space, see (pT]) and (0), as a factor. 

Once the amplitudes n < pqa \ U > are determined, the matrix elements in ([19| ) and 
(pip can be evaluated by quadratures in the manner described in H] and references therein. 



IV. THE TRANSITION OPERATOR 

On top of the vr-induced transition potential described in section || we include exchange 
potentials driven by 77, K, p, uj and K* mesons. The explicit expressions for the weak and 
strong Hamiltonians can be found in Ref. 0]. 
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The resulting one boson exchange expression in a nonrelativistic reduction for the pion 
is given in (^) which we rewrite here as: 



V^{q^ = -Gprni^ [ A + ^(Ti ■ qA (35) 

where g = (y^NNvr is the strong couphng constant for the NNvr vertex, /i is the pion mass, 
stands now for the momentum carried by the exchanged meson and the operators A and B 
contain the isospin dependence of the potential 

A = A^ Ti-T2 

B = B^n-r2. (36) 

For pseudoscalar mesons different from the pion we have an expression analogous to 
but making the following replacements: 



9 - 




- 




A - 




B - 


Br,. 



(37) 

when considering the exchange of the isoscalar 77-meson, and 

9 ~^ 9ank, 



A^l^ + Dl^ + ^n-r^]^ 

B^r^ + D-^ + ^n-rA (38) 



for the isodoublet kaon. 

The factor Mj^/M corrects for the fact that the nonrelativistic reduction of the strong 
ANK vertex is now proportional to giving a factor 1/M instead of 1/Mj 



N- 



In the case of vector mesons as the p, one obtains 0]: 



e{Fi + F2),^ ^\ 1 



with /i = mp, Fi = g^j^p, F2 = g'^^^ and where the operators a, (3 and e have the following 
structure: 

a = Up fi ■ T2 
P = Ppn-T2 

e = epTi-T2 (40) 
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We can get the nonrelativistic potential corresponding to the exchange of the rest of 
vector mesons by making the following replacements in Eq. (|39|): 



/i m^, 
Fi 

Fo 0^ 

^ ^ e^. (41) 

for the exchange of the isoscalar cj, and 



?ANK* ' 

F2 ^ - 



ANK* ' 

/^PC.V ^PC,V 

a + ^K- + ^—^1 ■ ^2, 

/^PC,T ^PC,T 

P ^ — ^ \- J->K* H ■ r2, 

(42) 

for the isodoublet K*-meson. 

In configuration space the potential including the exchange of all the mesons can be cast 
into the form: 

^(r) = E E = E E (r)oJ!f> 

i a i a 

= E \yl?ir)i^ + V^{r)a, ■ ajg + V? {r)S,,{r)i^h 

i 

+ (n'a^ . f + (1 - Pi x ■ r) V#(r)4t] , (43) 

where the index i runs over the different mesons exchanged {i = 1,...,6 meaning 
71, p,K,K* ,ri,Lj) and a over the different spin operators denoted by C (central spin inde- 
pendent), SS (central spin dependent), T (tensor) and PV (parity violating). In the above 
expression, particle 1 refers to the A and = 1(0) for pseudoscalar (vector) mesons. For 
isovector mesons (vr, p) the isospin factor ij^^ is ti ■ 72 for all values of a, for isoscalar mesons 
{ri,u) this factor is just 1, and for isodoublet mesons (K,K*) there are contributions pro- 
portional to 1 and to fi ■ T2 with coefficients that depend on the coupling constants and, 
therefore, on the spin structure piece of the potential denoted by a. 
For K-exchange the isospin factors in (^) are 

= 

/^PC /^PC 

ipv = ^ + D7 + ■ T2 (44) 
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and for K* exchange they are 



^PC,V ^ 
P(6) _ '-^K* V ^^ 



2 



PC,V 
PCV I "-^K* 



-Ti ■ T2 



f^PC,V I ^PC,T\ /^PC,V , ^PC,T\ 

^SS — — n r (^i^K* + ^k" ) H ^1 ' ^2 



2 

f(6) _ '-"k* I r-)PV I '-^K* 



(45) 



The different pieces with a = C, SS,T, PV, given in Ref. 0, are reproduced here 
for completeness 



-6{r) 



3 3 
1 + — + 



Attt 



1 + — =K';>yVpv{r,fi,) . 



(46) 
(47) 

(48) 

(49) 



where fii denotes the mass of the different mesons. The expressions for K^^ , which contain 
factors and couphng constants, are given in Table |. The explicit values of the strong and 
weak coupling constants, taken from Ref. [0], are shown in Table ||. 

Including monopole form factors Fj(g ^) = (A^ — fif)/{Af + q"^) at both vertices, where 
the value of the cut-off, Aj, depends on the meson (see Table |lD, leads to the following 
regularization of the potential: 



Vc{r; /i,) Vc{r; fi,) - Vc{r- A,) - A 



2 „-A,:r 



47r 



A,;r 



Vssir; ^^^) Vss{r; ^^^) - Vssir; Ai) - A.. 



AC - f^i e 



2 „-Air 



Vrir; /x,) ^ VV(r; /x,) - Vrir; A,) - A. 
Vpvir; Hi) Vpv{r; Hi) - Vpv(r; A^) - 



A/ - Hi e 



An 

2 o-A,r 



A,T 



2 47r 
A,' - e-^^^ 
2 47r 



1 + 



A,r 



(50) 
(51) 
(52) 
(53) 



where Kf(r; Aj) has the same structure as V^(r; Hi), defined in Eqs. (|46|)-(^9D, but replacing 
the meson mass Hi by the corresponding cutoff mass Aj. 

The last step is the transition into the momentum space partial wave representation 
(1^). Of course we could have derived that directly from (p5 H39|) using the standard helicity 
formalism [ p^ . 

This leads to 



2 ./;/ 



TT 



<p{ls)j I \p{l's')j > = r~ drr^ji{pr)V^\r)jv{pr)5iv5ss' + 
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r drr^ji{pr)vi{r)ji,{p'r) < {Is)] | <?! • ^2 I il's')] > + 

TT JO 



2 



oo 



-i^'-'> / drr'ji{pr)V^'>{r)ji,{p'r) < {ls)j \ S^if) \ {l's')j > + 

71 Jo 



2 

TT 



drr 3i{pr)Vf!f{r) iv{p'r) < {ls)j | as ■ r | {l's')j > + 



2 (1 _ r c?rr2j,(pr)\/i(J(r) jHpV) < {ls)j \ [a^ x a^] ■ f \ {l's')j > (54) 



TT 



The radial integrations were carried out numerically. The angular momentum parts are 
standard and are given for instance in 0. 

Since our results show a strong dependence on the different meson contributions with 
varying signs, we would like to display the radial shapes of the four types of potentials 
(C, SS, T, PV) split into the different meson contributions. This is shown in Figs. ^ and 

for the central spin-independent, in Figs, ^a and for the central spin-dependent, in 
Figs. ^ and for the tensor and in Figs. ^ and §3 for the parity violating channels. Note 
that we have represented r'^V^^\ where V^^^ is the potential regularized by the effect of form 
factors, and that the expectation value of the isospin factor for each meson and channel has 
also been included. As expected, we observe that the vr-meson contribution is by far the one 
of longest range. More interesting is to note that, compared to the pion, all the other mesons 
play a relevant role in a wide range which extends up to about 1.5 fm. On the right hand 
side of the figures we have plotted the full potential obtained when all meson contributions 
are added. We observe that the full potential is clearly different from the vr-only one, shown 
on the left hand side of the figures. We also see that in the spin independent central channel 
(Figs. ^,b) the contribution of the vector mesons compensate each other giving rise to a 
practically negligible transition potential for both isospin channels. 

In Figs. |^,b only the vector mesons appear since they are the ones that contribute to 
the spin-independent channel. These figures show that the three potentials have about the 
same range and their contribution is similarly relevant. 

As seen in Figs. ^,b, the K* meson gives a very important contribution to the central 
spin-dependent channel. We also observe that, except in the intermediate range where the 
potentials change sign, there is a constructive interference between the pseudoscalar and 
vector components of each isospin-like pair [ (7r,p) , (K,K*) , {ri,Lj) ]. Note that, in the 
T = 1 channel, the tu-meson lies very close to the p-meson potential. This is due to the 
similar value of the p and u masses and to the fact that the combination of strong and weak 
coupling constants building up i^"ss^ (see Tables | and |I| ) gives, by chance, a very similar 
value. This behavior is not observed for the T = channel because, due to its isovector 
character, the p-meson contains an additional factor of —3 compared to the ci;-meson, as can 
be clearly seen in Fig.^a. 

The tensor transition potential is shown in Figs |^a,b. In this case, we observe a de- 
structive interference pattern for each pair of isospin-like mesons. In the T = 1 channel the 
K*-meson clearly stands out with respect to the other mesons and, for the same reasons 
explained above, the p and u contributions are again very similar. 
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Figs. |^,b show the parity violating contributions. Here, the interference is constructive 
for the {r),u}) pair and destructive for the (7r,p) pair in both isospin channels. The (K,K*) 
pair shows a destructive interference in the T = channel and a constructive one in the 
T = 1 channel. In these figures the longest range of the pion contribution stands out quite 
clearly over the other mesons, especially in the T = channel. 



V. RESULTS 



We used a hypertriton wave function based on the Nijmegen 93 NN potential |jT3| and the 



Nijmegen YN interaction 0], which include the A — S transitions. The number of channels 
(see section pT|) used in the solution of the corresponding Faddeev equation is 102. This 
leads to a fully converged state, which has the proper antisymmetrisation among the two 
nucleons built in. Also the NN and YN correlations are exactly included as generated by 
the various baryon-baryon forces (see The UNN part of the state has a probability of 
0.5 % and will be neglected. 

The deuteron and the final state interaction among the three nucleons is generated using 
the Nijmegen 93 NN force, including the NN force components up to total two-body angular 
momentum j = 2. This is sufficient to get a converged result for the nuclear matrix element. 

Since the total three-body angular momentum is conserved, the Faddeev equation (|20| ) for 
the final state interaction (FSI) has to be solved only for total three-body angular momentum 
J = |, but for both parities due to the parity violating transition potential. 

The total nonmesonic decay rate turns out to be F = 0.21 x 10^-, which is 0.55 % of the 
free A decay rate. Fa = 3.8 x 10^^. This is about one order of magnitude smaller than what 
has been found in the rough estimate [^, which was based only on vr-exchange, a simplified 
hypertriton wave function and the absence of FSI. In the following, we show that the final 
value for the total decay rate comes from many dynamical ingredients, which all contribute 
significantly. Therefore, that quantity will be an important test for our understanding of 
that system and should be measured. 



Table |I| shows the individual contributions of the six mesons to Tnm and the way each 
meson contributes to the final result. We see that the pion by itself provides the largest 
contribution, followed by K, K* and u. Adding the meson contributions one by one can 
yield a strongly varying sequence as seen in Table ^ choosing a special but arbitrary order. 
In view of Figs. ^|-^ this is hardly surprising. The final result, however, is such that one 
ends up close to the value with pion exchange only. 

The total decay rate is the sum of the partial rates for the nd and nnp decays. Our 
result for the pion only are shown in Table |V[ The nnp contribution is clearly dominant. 
We also show the plane wave impulse approximation (symmetrized) (PWIAS) results. They 
are defined by evaluating the nuclear matrix elements in Eqs. ([T9|) and (^) keeping only the 
first terms, respectively. The comparison of PWIAS to the full result (keeping both terms in 
Eqs. ( [T9|) and (^If )) underlines the importance of the final state interaction, which reduces 
the rate. Finally the parity conserving and parity violating contributions are listed and it is 
seen that they are comparable to each other, though with a slight dominance of the parity 
conserving part. 

The corresponding numbers including all mesons are also displayed in Table Again 
the final state interaction is very important and reduces the PWIAS results by about a 



14 



factor of 2. Now for all mesons included the parity conserving part is clearly dominant. 

There is often a separation of p- and n-induced decay in the literature [§,0. They act 
clearly coherently and strictly spoken cannot be separated experimentally. Theoretically, 
however, we can choose in the intermediate state in front of | \l/3j:^ > in Eqs. ([19|) , (|20|) and 



(pl|) a situation that the A-particle chooses either a proton or a neutron as its partner for 
meson exchanges. This amounts to keeping only the first or second term on the right hand 
side of Eq. (^), respectively. As already mentioned above this requires to keep also T = 3/2 
states in the final state. 

Our results for the pion only and for all 6 mesons are displayed in Table |V]. For the nd 
breakup clearly the separate n- and p-induced decay rates do not add up to the total nd- 
decay rate, which tells that there is interference. On the other hand for the nnp breakup the 
separate contributions from the n- and p-induced decays add up to the total nnp-decay rate. 
However, this does not imply automatically that they can be separated experimentally. We 
come back to that interesting issue below. We see that for nd and nnp decay the p-induced 
decay is stronger. 

Focusing on the PWIAS for the nnp breakup of Eq. (|I9D one has three contributions 



< (t)l I (1 + P)0(1,2) \^iH>= -^il< pmim2ViV2 \ 3< qm^v^ \ 0(1,2) | > 



+ 23 < prnim2 z/iZ/2 | i < qrrisUs \ 6(1, 2) | ^3 j:^ > 



+ 31 < pmim2 z/iZ/2 I 2 < qm^i'^ \ 6(1, 2) | ^!z^h > (55) 

In the above equation we applied the P operator to the left. As always, the subscripts refer 
to the particles in states with momenta p, q, spin magnetic quantum numbers m and neutron 
or proton labels u. Using Eq. ( pO] ) the isospin matrix elements can simply be calculated with 
the result that in the first matrix element on the right hand side of Eq. (^5]) nucleons 1 and 
2 are two neutrons for n-induced decay and a neutron-proton pair for the p-induced decay. 
Also that first matrix element peaks at g*= 0, which means that nucleons 1 and 2 share the 
total energy and emerge back to back. Under this kinematical condition the other two matrix 
elements are strongly suppressed, as is manifest if one expresses the momenta occurring in 
these two matrix elements in terms of p and g* = of the first matrix element. If we denote 
the p from the first matrix element as pi2, then it turns out that p = — |pi2(— |pi2) and 
q = Pui—Pu) in the second (third) matrix element, respectively, and for such a g- value 
I ^3 j:^ > is suppressed. 

The other two matrix elements also peak if particles 1 or 2 emerge with zero momenta. 
Therefore, we have to expect 3 peaks. Let us now take a closer look at the quantity 



dr^'^'^'^P / dpd6 , defined in Eq. |T0|. That quantity, suitably restricted to certain or all meson 
exchanges, summed over all p and 6 values provided the various values of Tables fT|-|V[ For 
the choice of Jacobi momenta (^) the 3 peaks are located as sketched in Fig. |^. Energy and 
momentum conservation requires that + = p^ax' where Pmax is the maximal p-value. 
For the available energy Pmax ~ 2 fm~^. As an example we regard ^2 = 0. Then p = ^ki and 
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q = —ki. It follows that ki = p^ax and consequently p = ^Pmax- For ^3 = the momenta ki 
and ^2 have to be back to back and there will be no 9 dependence; for /c2 = the momenta 
^3 and ki are opposite to each other, therefore 6 = tt and finally for ki = the momenta 
/c3 = —k2 and ^ = 0. For PWIAS the result is displayed in Fig. |^ for the exchange of all 
mesons. The corresponding results for vr exchange only is qualitatively similar but larger by 
about 50 %. As expected we see the 3 peaks at the proper locations. The variation with 
6 for the maximal value is due to the factor sin^ in the expression ([T0|). Note that the 
sin 6'-dependence also removes the highest peak values for all 3 peaks. 

For the choice of nucleons 1 and 3 being neutrons and 2 a proton we have thus to expect 
that for a neutron induced decay there should be only one peak at the position p = |pmax 
and 6 = n, which is indeed the case as shown in Fig. ^. Note that in evaluating the nuclear 
matrix element of Eq. |10| we fixed the isospin magnetic quantum numbers z/i, z/2, z/3 to 
be — |, |, — |. This refers to all figures 3^l3|. For p- induced decay we expect two peaks 
corresponding to either fci = or ^2 = 0. And this is what comes out and what is shown in 
Fig. |TU[ Regarding Figs. P-pHj we see that the areas populated by n- and p-induced decays 
appear to be well separated in phase space and seem to add up essentially incoherently to 
the full result. A closer inspection, however, will be carried through below, which leads to 
a different result. 

Now let us turn on the final state interaction. For 3 nucleons interacting among each 
other one knows from Nd breakup reactions 0, that cross sections are strongly enhanced 
if two nucleons emerge with equal momenta. This is due to the strong interaction in the 
""^^o state, where the NN t-matrix has a pole close to zero energy (virtual state). These 
enhancements will be called final state interaction peaks (FSIP) in the following. Fig. |^ 
shows the positions, where this happens in the p — 9 plane. For the case p = clearly no 
6'-dependence is present. For all meson exchanges Fig. |ll] shows dV'^^'^^^ / dpdO including 
the full final state interaction. We see again the 3 peaks already known from the PWIAS 
result, but with reduced heights according to the already known reduction of the rate due 
to the final state interaction. In addition there are two more little peaks caused by the final 
state interaction for two pairs of nucleons, where the nucleons forming a pair have equal 
momenta, respectively. The final state interaction peak for the third pair is suppressed by 
the kinematical factor p^ in (|TU]). The p- and n-induced pictures (Figs. |T^]-[T^) keeping 
the full final state interaction again look qualitatively similar to the ones evaluated in the 
PWIAS approximation, only the final state interaction peaks are added. 

It is interesting to see despite the fact that FSI decreases the dV"^^^^"^ / dpdO values sig- 
nificantly, which means a strong rescattering among the three nucleons, that there is only 
one peak for the n-induced decay and the rescattering does not populate the other two peak 
areas. Corresponding is also true for the p-induced decay. Again the events for n- and 
p-induced decays seem to add up incoherently in the quantity dV'^^'^^'^ / dpdO. 

Let us now discuss the form (|lTD of the differential decay rate expressed in individual 
momenta in the total momentum zero frame of the decaying hypertriton. After averaging 
over the initial state polarization and summing over the spin magnetic quantum numbers of 
the final 3 nucleons the decay rate can depend only on the angle between the two nucleon 
detectors, 9i2. We now show two sets of figures. Figs. p!^-[15|, for various Oi2's 0° < ©12 < 
180° and for pion exchange only. In Fig. we compare PWIAS to the full calculation 
including final state interaction. Thereby the two detected nucleons can be either a proton 
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neutron pair or two neutrons. For 0i2 = 180" we see a strong enhancement for S ~ 110 
MeV. (This corresponds to a location around the middle of the locus in Fig. |^.) In PWIAS 
this is caused by the fact that there the wavefunction enters at g = /C3 = 0. Like in 
Nd scattering we shall from now on call such a configuration, where one final nucleon has 
zero momentum, a quasi free scattering (QFS) configuration. The final state interaction 
reduces that enhancement, but it is still pronounced. In addition we see two FSIP's in 
the full calculation. (They have to be absent of course in PWIAS). For B12 = 160° that 
enhancement is reduced and two peaks emerge at the beginning and the end of the S- 
curve. Since there either Ei or E2 are small, we have again configurations, which are close 
to QFS conditions, now for the nucleon pairs 2,3 and 1,3, respectively. This explains the 
additional enhancements. Now at G12 = 120° the enhancement in the middle of the S-curve 
has disappeared. That point on the S-curve corresponds exactly to the so called space-star 
configuration in aN + d^N + N-\-N process. All three nucleons receive the same energy 
and emerge completely symmetrically under 120° pairwise angles. This is far away from 
QFS conditions and no enhancement is expected. Of course, at the beginning and end of 
the S-curve the peaks correspond again to conditions close to QFS. The situation remains 
similar at 0i2 = 90° and 60°. Finally a new structure appears at 612 = 20° and above all at 
©12 = 0° in the middle of the S-curve. This is a FSIP, which is fully developed for 612 = 0°. 
In principle it could be used to extract information on the np and nn scattering lengths like 
in Nd breakup processes. 



In Fig. 15 we compare the full calculations to the separate decay rates for n- and p- 



induced processes. Note that in this figure PN means that the proton is nucleon 1 and a 
neutron is nucleon 2, which corresponds to z^i = |, 1^2 = ~| and 1/3 = — | in the matrix 
element of Eq. (|Tl|) . 

At G12 = 180° the decay rates in the center of the S-curve are essentially given by the 
p-induced process if a proton neutron pair is registered and by the n- induced process if 
two neutrons are registered. Already at 612 = 160° this is no longer true. The p-induced 
rate and even more the n-induced rate is significantly lower than the rate built up by the 
full physical process. It is interesting to see in cases of Fig. |l^(b-f) that the n-induced 
decay rate for pn detection at the upper end of S is practically identical to the decay rate 
corresponding to the full physical process. The reason is that at the upper end of 5* the 
energy Ei (the proton energy) is nearly zero (see Fig. Q), and therefore two neutrons carry 
essentially all the energy. Without FSI such a case can only be generated by a n-induced 
decay and quahtatively this picture does not change due to FSI. In case of p-induced decay 
its rate for pn and nn detection at the lower end of S is practically identical with the full 
decay rate. The reason is similar as for the n-induced decay. At the lower end of S the 
energy E2 (a neutron energy) is nearly zero. Therefore a proton-neutron pair has to carry 
essentially all the energy and this has to be generated by a p-induced process. All these 
enhancements at the lower and upper end of S are QFS-like cases. 

Interesting are also the FSIP's, especially pronounced at 612 = 0° and 180°. At 612 = 0° 
and for neutron-neutron detections the proton has to fly in the opposite direction, therefore 
the p-induced process has to be mainly responsible for the peak, as is the case. For p-n 
detection, however, both p- and n-induced decays can contribute to a FSIP and they do. 
Apparently the p- and n-induced decays have to interfere, since the individual rates do not 
add up to the total physical decay. As a further example we comment on the left strong FSIP 
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at 012 = 180°. There are two peaks according to pn and nn detections. For nn detection (the 
neutrons have opposite momenta) this has to be necessarily a pn FSIP. As we see from the 
figure it receives contributions from n- and p-induced processes, again coherently. In case 
of pn detection (the proton an a neutron have opposite momenta) the special location on 
the S'-curve (high proton energy) requires that it is a nn FSIP. Therefore the very dominant 



contribution has to come from the p-induced process. This is clearly visible in Fig. 15 



The corresponding curves, when all mesons are included, are qualitatively the same and 



only very few examples are displayed in Fig. |1^. The heights of the FSIP's have changed, 
however, significantly. 

Finally we consider the question, whether the total n- and p-induced decay rates can be 
separated experimentally. As we already saw this is not possible in the nd decay channel. 
There the two processes interfere coherently and the individual theoretical rates do not sum 
up to the total rate. In the 3N decay channel the total rate is very close to the sum of the 
individual rates for the n- and p-induced processes. Also in the p — 6 representation of 



dT (see Eq. (Rw) and displayed in Figs. 12-13 the events from the two different processes 



appear to be nicely separated. On the other hand, in the 0i2 — S representation, which 
is directly accessible using two detectors, we saw cases where an interference was manifest. 
The total decay rate into 3 nucleons is 

f d^Cr ^ o /■'^ /■Smax(012) 

F'^+'^+P = / ^^^^ dhdk2dS = Srr^ rfGiasineis/ ci5a(0i2,5) (56) 

J dkidk2dS Jo Jo 

We used the fact that d^a/dkidk2dS depends only on 9i2 and S and introduced the length 
'S'max(0i2) of the S-curve depending on 612. First of all we notice immediately that the pure 
QFS cases, where a final nucleon momentum is zero, do not contribute, since for ^3 = 
012 = 180° and for fci = or ^2 = the phase space factor in Eq. (|TT]) is zero. Nevertheless 
an angular configuration with 0i2 = 180° can and should be measured by itself, since along 
the S-curve there will be one point with the exact QFS conditions and as we saw in Fig. p!5 
there the n- and p-induced processes can be cleanly separated. One measures either a pn 
or a nn pair and they are generated by p- and n-induced decays, respectively. 
We now discuss the quantity 

/■<S'max(012) 

7(012) = 87r2 sin 012 / d^a(0i2,5) (57) 

Jo 

for np and nn detection, respectively. This is shown in Figs. p!7|-p!^ together with the individ- 
ual contributions of the p- and n-induced processes. We see strong peaks near 0i2 = 170° 
for the full processes. The corresponding values for the p(n)-induced decay are similar in 
the peak area for pn(nn) detection, while the n(p)-induced values are small. 

At smaller angles 0i2 the p- and n-induced values are similar to each other in case of pn 
detections, while for nn detection the p-induced quantities dominate. A closer inspection 
reveals that the sum of the p- and n-induced values for each 0i2 do not sum up very well 
to the value according to the true physical process, but there are differences up to 10 % (nn 
detection) around 0i2 = 170°. This is a clear signal for interference. 

Let us quantify this question. The representation (p!0| ) of dV"'^'^^^ / dpdO , which has been 
displayed in Figs. p|-[T3| yields the decay rates for the individual n- and p-induced processes 
as well as for the full physical process (neutron and proton induced) when integrated over 
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the whole p-6 plane. According to Figs. [T^ and |T3| the n- and p-induced processes appear 
to receive their contributions from well separated areas in the p-9 plane. Quantitatively, 
however, this is not true. Restricting the integration in 9 and p to the region where the peaks 
in Figs. p!2|-|r^ are located results in only a certain fraction of the full rates. Quantitatively, 
if we fix that fraction to 60 %, say, for p- and n-induced decays, respectively, we find that 
the regions displayed in Fig. |T9| contribute. In choosing a certain fraction we always start 
integration from the highest values (located in the peaks) downwards and stop when the 
assumed fraction has been reached. Except for a small domain (6* ^ vr and p large) the 
regions for p- and n-induced decays are clearly separated. Of course for fractions smaller 
than 60 % this will be even more the case. For fractions larger than 60 %, however, the 
regions overlap considerably. An example for 90 % is also shown in Fig. 0. Clearly in such 
a case the events coming from for p- and n-induced decays cannot be separated any more 
exp er iment ally. 

We show the fractional decay rates evaluated over corresponding increasing regions in 



Tables 0, [VI I| and in Figs. p!^-pO|, where the results refer to full calculations and include 
all meson exchanges. There the p- and n-induced rates are compared to the observed rate 
Fphysicai, produced by the full physical process. Thus if we require that the measured value is 
equal to the p- or n-induced decay within a few percent, one has to restrict the integration 
in the 9-p plane to certain subregions and the rates to only about 60 % of the full rate. 

Let us map the 9-p values into the variables 012-5*, which are directly accessible ex- 
perimentally. This is shown in Fig. ^ for the 60 % and 90 %. That picture refers to the 
detection of a neutron (particle 1) - proton (particle 2) pair. Fig. ^ tells that n-induced de- 
cay can be found under all O12 angles for small S'-values. {E2 is then small.) The p-induced 
decay on the other hand is to be found for all the Oi2-angles around maximal S'-values and 
in the region 160° < O12 < 180° for medium-large S'-values. (For small neutron energies.) A 
correspondingly modified figure could be shown if two neutrons would be detected. If on the 
other hand 90 % of the corresponding rates should be detected then the two detectors would 
receive events from both processes under the same angle O12 and for the same energies in 
a large portion of the phase space. Thus experimentally it is not possible to separate those 
processes. 

We have to conclude that despite the fact that r"'"*"""'"^ is rather close to the sum of 
r"p^"^^ and r"^"'"''^, the latter values cannot be determined experimentally, only fractions, 
the smaller, the cleaner. 



VI. SUMMARY AND CONCLUSIONS 

For the first time the nonmesonic hypertriton decay has been calculated based on rigorous 
solutions of 3-body Faddeev equations for the hypertriton and the 3N scattering states of 
the final 3 nucleons. Realistic NN and hyperon-nucleon interactions have been used. In the 
meson exchange process the pion exchange is dominant, but the other included mesons 77, 
K, a;, p and K* provide also significant contributions of various signs and magnitudes. The 
final state interaction turned out to be very important and reduces the rates for PWIAS 
by about a factor of 2. The total nonmesonic decay rate turns out to be 0.55% of the 
free A decay rate and is smaller by an order of magnitude than a previous estimation [4] 
which used the pion-exchange model, a much more simplified hypertriton wavefunction 
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and no FSI. While the p- and n-induced decays add up manifestly in a coherent manner 
in the nd decay process, the nnp decay rate is rather well given as the sum of the n- 
and p-induced decay rates. Nevertheless these individual decay rates cannot be measured 
separately. Only fractions thereof can be obtained, when the contributions arise from non- 
overlapping regions in phase space. This subject has been thoroughly discussed in section |V]. 
Detailed information has been given regarding the separation of p- and n-induced differential 
decay rates experimentally. 

Since the decay rates depend sensitively on the number and type of mesons exchanged, 
it will be an interesting testground for the dynamics of these meson exchanges, which are 
driven by weak and a strong vertices. At the same time the decay rates probe the hypertriton 
wavefunction and the reaction mechanism of the three outgoing nucleons through their strong 
final state interaction. The latter one is especially manifest in the strong FSIP's, where two 
nucleons leave with equal momenta. 

The evaluation of the pionic decay into the various bound states and continuum channels 
can be performed in a similar manner and is planned. 
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FIGURES 



FIG. 1. Locus for kinematically allowed events in the E1-E2 plane and Q12 = 180° together 
with our definition of the choice for the arclength S = 0. From that point on S is evaluated for 
each point on the locus in the counterclockwise sense. 



FIG. 2. The vr-induced AN NN exchange process of Eqs. (23-24) 



FIG. 3. Individual meson contributions to the isoscalar (a) and isovector (b) central spin inde- 
pendent regularized potential r^V^(r). The potential obtained by adding all meson contributions 
for the isoscalar (c) and isovector (d) case. 

FIG. 4. Individual meson contributions to the isoscalar (a) and isovector (b) central spin de- 
pendent regularized potential 'f'^Vgg{r). The potential obtained by adding all meson contributions 
for the isoscalar (c) and isovector (d) case. 

FIG. 5. Individual meson contributions to the isoscalar (a) and isovector (b) tensor regularized 
potential r^Vjf(r). The potential obtained by adding all meson contributions for the isoscalar (c) 
and isovector (d) case. 

FIG. 6. Individual meson contributions to the isoscalar (a) and isovector (b) parity violating 
regularized potential r^Vpy(r). The potential obtained by adding all meson contributions for the 
isoscalar (c) and isovector (d) case. 



FIG. 7. The location of the 3 peaks corresponding to ki = (i=l,2,3) and the FSI peaks in the 
9-p plane (see text). 

FIG. 8. The differential decay rate dT^~^^^^ / dpdO for PWIAS and exchange of all mesons. 

FIG. 9. The n~induced differential decay rate dT^^~^~^^ / dpdO for PWIAS and 7r-exchange only. 
Nucleons 1, 2 and 3 are chosen to be a neutron, a proton and a neutron. 

FIG. 10. The same as in Fig. ^ for p-induced decay. 

FIG. 11. The differential decay rate dT'^'^'^'^P / dpd9 with full inclusion of the final state inter- 
action and exchange of all mesons. 

FIG. 12. The n-induced differential decay rate dV^^-^^^ / dpdO with full inclusion of the final 
state interaction and vr-exchange only. 
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FIG. 13. The same as in Fig. 
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for p~induced decay. 



FIG. 14. The differential decay rate dT'^'^'^'^P / dQi2dS for various angles 0i2 and vr-exchange 
only. The two detected nucleons are either a p(particle l)n(particle 2) pair or two neutrons. PWIAS 
is compared to the treatment including the full final state interaction. 

FIG. 15. The n-and p-induced differential decay rates dF^^-^'^^^-l' /dQudS in comparison to the 
physical rate for various angles 0i2 and vr-exchange only. The two detected nucleons are either a 
p(particle l)n(particle 2) pair or two neutrons. The final state interaction is fully included. 



FIG. 16. The same as in Fig. 15 for the exchange of all mesons. 



FIG. 17. Differential decay rates integrated over the S-curve as a function of ©12. The p- and 
n-induced cases are compared to the physical process. The final state interaction is fully included 
and all mesons are exchanged. The curves belong to the case that a proton and a neutron are 
detected. 



FIG. 18. The same as in Fig. 17 for the case that two neutrons are detected. 



FIG. 19. The separate regions in the 9-p plane contributing to (a) 60 % and (b) 90 % of the 
rates of n- and p-induced decays. Note the strong overlap of the different processes in phase space 
in case (b). 

FIG. 20. The separate regions in the ©i2-»S' plane contributing to (a) 60 % and (b) 90 % of 
the rates of n- and p-induced decays. Note the strong overlap of the different processes in phase 
space in case (b). (Particle 1 is a neutron and particle 2 is a proton.) 
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TABLES 



TABLE L Constants appearing in the weak transition potential for different mesons (in units 
of Girm^^). The strong and weak couphng constants have been taken from Ref. 





K'ss 






vr 


TT y IN IN TT 


TT y IN iN TT 


^ y IN IN TT 

"2Mjv 


2M 2Mn 


2M 2MAr 


r] 


-Biyj 0]\n\rr) 

2M 2Mn 


// IN INI/ 

2M 2Mjv 


y* L\ l\ 1 j 

''2M^ 


K 


1 Sank 
2Mjv 2M 


-1 

J- 5ank 
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P Onnp'^P 


2M 2M7V 


Q^p + /5p 5nnp "1" 5nnp 

2M 2MAr 


Snnp ~I~ Snnp 

2M7V 


^ 9nnuj 




Q^a; ~l~ Pu) 9NNiu ~^ 9nnuj 
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2M 2Mn 


2M 2Mn 


" 2MAr 




o ^ Sank* "1" 5ank* 


1 5ank* + Sank* 


Sank* "1" Sank* 


2Mjv 2M 


2Mn 2M 


2Mn 
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TABLE II. Strong (Nijmegen) and weak coupling constants and cutoff parameters for the 
different mesons. The weak couphngs are in units of G^mTr^ = 2.21 x 10"''. 



Meson 


Strong c.c. 


Weak c.c. 
PC PV 


(GeV) 


TT 


9NNir — 13.3 
gA^n = 12.0 


B^=-7.15 


^^=1.05 


1.30 


V 


5NNr) = 6.40 
gAArj = -6.56 


5^=-14.3 


^^=1.80 


1.30 


K 


Sank = -14.1 
Snek = 4.28 


C^'==-18.9 
£»P«=6.63 


C^^=0.76 
£»P^=2.09 


1.20 


P 


9nnp = 3.16 
9nnp — 13.3 
9l^p = 
5Lp = 11-2 


Q;p=-3.50 


ep=1.09 


1.40 




„V _ -iri ir 

9nnui — ■'-U-'J 
Snno; — 3.22 
91ao. = 7.11 
ffL, = -4.04 


au)=—S.69 


6^^= —1.33 




K* 


sInk* = -5.47 
9Lk* = -11-9 
5nsk* = -3.16 
=6.00 


C^?'^=-3.61 
C^F=-17.9 

DP5-^=_4.89 

i:),''?'^=9.30 

Iv 


cpy=-4.48 

DPy=0.60 


2.20 



TABLE III. Decay rates in units of s for individual meson exchanges and for partially 
summed up exchanges. 



■pmeson 


■p partially summed 


TT 


0.2412 xlO^ 


TT 


0.2412 xlO^ 


7] 


0.4826 xlO^ 


TT + 77 


0.2299 xlO^ 


K 


0.5422 xlO^ 


TT + r] + K 


0.9267 xlO^ 


P 


0.7647 xlO^ 


TT + r] + K + p 


0.7502 xIQS 


UJ 


0.4372 xlO^ 


TT + r] + K + p + u 


0.1752 xlO^ 


K* 


0.5569 xlO^ 


TT + rj + K + p + u + K* 


0.2126 xlO^ 
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TABLE IV. Selected decay rates in units of s ^ for 7r-exchange only and for exchange of all 
mesons. 





TT-exchange only 


Exchange of all mesons 


p n+d 
^ PWIAS 


0.59 xlO^ 


0.47 xlO^ 


p n+d 


0.15 xlO^ 


0.22 xlO^ 


p n+n+p 
PWIAS 


n 4fi V 1 


n 'if, V 1 
u.ou A lyj 


■p n+n+p 


0.23 xlO^ 


0.19 xlO^ 


p n+d 
^ PC 


0.88 xlO^ 


0.22 xlO^ 


p n+d 
J- PV 


0.59 xlO^ 


0.22 xlO^ 


p n+p+p 
^ PC 


0.13 xlO^ 


0.12 xlO^ 


p n+p+p 
^ PV 


0.92 xlO^ 


0.73 xlO^ 



TABLE V. Proton- and neutron- induced decay rates in units of s for 7r-exchange only and 
for the exchange of all mesons in comparison to the total nd and nnp rates. 





TT-exchange only 


Exchange of all mesons 


pn+d 
^(") 


0.20 xl{f' 


0.64 xlO*^ 


pn+d 
^P) 


0.78 xlO^ 


0.70 xlO^ 


pn+d 


0.15 xlO^ 


0.22 xlO^ 


p n+n+p 


0.59 xlO^ 


0.57 xlO^ 


p n+n+p 
^(P) 


0.18 xlO^ 


0.13 xlO^ 


p n+n+p 


0.23 xlO^ 


0.19 xlO^ 
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TABLE VI. Fractional proton- and neutron-induced decay rates in units of s"-^ 
r(p),(n) = I '^^^p)"(n) integrated over subdomains in the 6-p plane, where proton induced decay 
dominates. They are compared to the corresponding fractional physical rate. 



area (P) 


^(p) 


r(n) 


^physical 


■^(p) /-^physical 


60 % 


0.80 xlO^ 


0.20 xlO^ 


0.82 xlO'^ 


98 % 


70 % 


0.94 xlO'^ 


0.48 xlO^ 


0.98 xlO^ 


96 % 


80 % 


0.11 xlO^ 


0.99 xlO^ 


0.12 xlO'^ 


92 % 


90 % 


0.12 xlO^ 


0.23 xlO^ 


0.14 xlO^ 


86 % 



TABLE VIL Fractional proton- and neutron-induced decay rates r(p) = / dF^"^"^^^ in units 
of integrated over subdomains in the 9-p plane, where neutron induced decay dominates. They 
are compared to the corresponding fractional physical rate. 



area (N) 


^(p) 


r(n) 


^physical 


■^(n) /-f^physical 


60 % 


0.34 xlO^ 


0.34 xlO^ 


0.35 xlO'^ 


97 % 


70 % 


0.11 xlO^ 


0.40 xlO^ 


0.47 xlO^ 


85 % 


80 % 


0.29 xlO^ 


0.45 xlO^ 


0.72 xlO^ 


63 % 


90 % 


0.50 xlO^ 


0.51 xlO^ 


0.10 xlO^ 


51 % 
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